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Abstract. We study transport properties of two Josephson junctions coupled by an 
external shunt resistance. One of the junction (say, the first) is driven by an unbiased ac 
current consisting of two harmonics. The device can rectify the ac current yielding a dc 
voltage across the first junction. For some values of coupling strength, controlled by an 
external shunt resistance, a dc voltage across the second junction can be generated. By 
variation of system parameters like the relative phase or frequency of two harmonics, 
one can conveniently manipulate both voltages with high efficiency, e.g., changing the 
dc voltages across the first and second junctions from positive to negative values and 
vice versa. 
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1. Introduction 

Transport phenomena in periodic structures by harvesting the unbiased external time- 
periodic stimuli and thermal equilibrium fluctuations have been one of the hottest topic 
in nowadays science. Examples range from biology and biophysics [lj to explain directed 
motion of biological motors [2] and particle transport in ion channels [3], or to design 
new separation techniques [4j, to meso- and nano-physics covering newest and up-to- 
date experiments with optical lattices [5], persistent currents in quantum rings [6] and 
Josephson junctions [TJ, to mention only a few. The physics of the latter has been 
studied for almost five decades now [SJ. Devices with Josephson junctions constitute 
a paradigm of nonlinear systems exhibiting many interesting phenomena in classical 
and quantum regimes. Yet new and interesting phenomena arise as researchers are able 
to use an advantage of powerful computer simulations followed by real experiments on 
junctions proving theoretically predicted findings j£l [101 ITT] . In this paper we study a 
system of two coupled Josephson junctions driven by an external ac current having two 
harmonics. Such driving has been considered as a source of energy pumped into the 
system and as an agent transferring the system to non-equilibrium states. Transport 
phenomena in non-equilibrium states are of prominent interest from both theoretical 
and experimental point of view. For example, transport processes driven by biharmonic 
driving have been intensively studied in variety of systems and in various regimes. We 
can quote: Hamiltonian systems |13j, systems in the overdamped regime [T2l IT4"] and 
for moderate damping [JJ)J [TB], dissipative quantum systems [T7], electronic quantum 
pumps [18J and quantum control of exciton states [19], cold atoms in the optical lattices 
[HI 121] , soliton physics [22] and driven Josephson junctions [23]- Here we intend 
to consider another aspect of transport in the system which consists of two coupled 
elements (subsystems). In particular, we take into consideration two coupled Josephson 
junctions. In a more general context it is a system which can be described by two 
degrees of freedom (it can consist of two interacting particles or one particle moving on 
a two-dimensional substrate) [24J. The question is: What new properties of transport 
can be induced by interaction between two elements and under what conditions can 
transport be generated in the second element if the driving is applied only to the first 
element. 

The layout of the paper is as follows: In section 2, we detail the system of two 
Josephson junction coupled by an external shunt resistance. Next, in section 3, we 
specify external driving applied to one of the junctions. In section 4, we analyze generic 
properties of dc voltages across the first and second junctions. In section 5, we elaborate 
on the voltage control. Section 6 provides summary and conclusions. 

2. Model and its dynamics 

Let us consider a system of two coupled resistively shunted Josephson junctions. The 
equivalent circuit representation of the system is illustrated in figure[TJ The junctions are 
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Figure 1. The system of two resistively shunted Josephson junctions coupled by an 
external shunt resistance R3 and driven by the currents I\(t) and -^(i)- 



characterized by the critical currents I c2 ), resistances (Rx,R 2 ) and phases (0i,0 2 ) 
(more precisely, the phase differences of the Cooper pair wave functions across the 
junctions). The junctions are externally shunted by the resistance R% and driven by 
external currents I\(t) and 12(f) ■ We also include into our considerations Johnson 
thermal noise sources an d £, 3 (t) associated with the corresponding resistances 

Ri, R2 and R3. 

From the Kirchhoff current and voltage laws, and two Josephson relations one can 
obtain the following equations 

R h 

= (R 2 + R 3 )[h(t) - I cl sin^i] + R 2 [h(t) - / c2 sin0 2 ] 
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- (R 2 + R 3 )Ut) - Ra&(t) - RsUt), 
= (Rx + R 3 )[h(t) - / c2 sin0 2 ] + #i[/i(t) 

- (Rx + R 3 ) 6(t) - Ri + Rs &(*), 

where the dot denotes a derivative of 0j = <pi(t) (i = 1,2) with respect to time t and 
R — Rx + R 2 + R 3 - We assume that all resistors are at the same temperature T and 
the thermal equilibrium noise sources are represented by zero-mean Gaussian white 
noises d(t) (i = 1,2,3) which are delta-correlated, i.e. (£i(t)£j(s)) oc SijS(t — s) for 
i,j G {1,2,3}. Note a symmetry property of equations ([lj with respect to the change 
Rx <-> R 2 . 

The limitations of the model Q and its range of validity are discussed e.g. in 
Ref. |25j . In particular, we work within the semi-classical regime of small junctions for 
which a spatial dependence of characteristics can be neglected, a displacement current 
accompanied with the junction capacitance is sufficiently small and can be ignored, and 
photon-assisted tunneling phenomena do not contribute. 
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The dimensionless variables and parameters can be introduced in various way. Here, 
we follow Ref. |26| and define the dimensionless time r as: 



where 



v °- Ic r 1 + r 2 + r 3 > h -^r~ (3) 

are the the characteristic voltage and averaged critical current, respectively. The 
corresponding dimensionless form of equations Q reads 

4>x = Ji (t) - J c i sin 4> x + a[I 2 (r) - J c2 sin <p 2 ] + vD rji (r) , 



2 = a /3 [/ 2 (r ) - J c2 sin 2 ] + a[h(r) - I cl sin0i] + \Jaf3D 772 (r), (4) 

where now the dot denotes a derivative with respect to dimensionless time r and all 
dimensionless currents are in units of I c . E.g., l c \ — > I c i/I c . The parameter 

is the coupling constant between two junctions and can be controlled by the external 
resistance i? 3 . The parameter (3 = 1 + (R3/R1). The noises i]i(t) and i] 2 (t) are zero- 
mean 5-correlated Gaussian white noises, i.e. {'r]i(' r )Vj{ T ')) = %^( r — r ') f° r hJ ^ {1; 2} 
and are linear combinations of noises £,i{t),i = 1,2, 3. The dimensionless noise strength 
is D = AeksT /hl c , where ks is the Boltzmann constant. 

Note that the dimensionless form ^ is not symmetrical with respect to the change 
i?i -H- i?2- It is because of the definition of the dimensionless time ^ which is extracted 
from the first equation of the set ([!]) and, in consequence, the asymmetry of Vo in 
equation ^ with respect to Ri and R 2 . 

3. System driven by biharmonic current 

The dynamic behavior of two Josephson junctions described by equations Q is 
analogous to the overdamped dynamics of two interacting Brownian particles moving in 
spatially symmetric periodic structures and driven by external time-dependent forces. 
However, it is not a potential system. If we think of (pi and <p 2 as position coordinates of 
two particles, then <pi and <p 2 correspond to the velocity of the first and second particles, 
respectively. In this mechanical analogue, junction voltage is represented by particle 
velocity and time-dependent current Ij(r) is represented by external force. We can use 
this analogue to visualize junction dynamics. The main transport characteristics of such 
mechanical systems are: the long-time averaged velocity Vi = of the first particle 
and the long-time averaged velocity v 2 = (<fi 2 ) of the second particle. The question is 
under what circumstances one can induce transport characterized by directed motion 
of both particles in the stationary regime, i.e. when v 1 7^ and v 2 7^ 0. In terms of the 
Josephson junction system Q it corresponds to the dimensionless long-time averaged 
dc voltages v% = and v 2 = (<j) 2 ) across the first and second junction, respectively. 
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Transport in symmetric periodic systems like Q can be trivially induced by biased 
external drivings ii(r) or/and ^{t)- A non-trivial case is when zero-mean (non-biased) 
driving can induce transport. An example of zero-mean driving is the simplest harmonic 
signal /(r) = Acos(lot). From the symmetry property of equations Q it follows that 
such a current source is unable to induce non-zero averaged voltages, for details see 
[HI EH EH EH EH E9]. However, when the system is driven by external currents having 
two harmonics, i.e. when (3U] 

Ij(t) = Ai [cos(c^t) + eiCos(QiT + 0*)] , i = 1, 2, (6) 

then transport can be induced within tailored parameter regimes. 

In the paper, we will study a particular case when the biharmonic ac current ^ is 



0=0 e=n/2 



6 = 
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Figure 2. (color online). Illustration of the biharmonic current (|7| for a\ = 1 and 
selected values of other parameters. Only symmetric (8 — 0,tt) and antisymmetric 
(8 = tt/2, 37r/2) forms of the current are displayed for specific values of the relative 
phase 8. For clarity, the curves are offset vertically from one another to facilitate 
comparison. 
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applied only to one (say, the first) junction, and has the form 

h(r) = ai [cos(W) + e x cos(2wr + Q)\ , J 2 (r) = 0. (7) 

It is known that without coupling to the second junction (i.e. when a = 0), the 
averaged dc voltage v 2 across the second junction is zero while the dc voltage v i across 
the first junction is non-zero in some regions of the parameters space. We can formulate 
several fundamental questions: (a) What is the influence of coupling a on the voltage 
v i? (b) Can the voltage across the second junction be induced by coupling ^ and how 
it depends on its strength? (c) Can the voltages reversal be obtained by changing the 
control parameters? (d) Can the voltages v\ and v 2 have opposite signs? Properties 
of the ac current I\{t) determine whether the voltages V\ and v 2 display the above 
mentioned features. One can distinguish two special cases of the two-harmonics ac 
current I\{t): 

(i) The first case is when there is such To that I\(tq + t) = I\(tq — t). It means that 
the driving is symmetric or invariant under the time-inversion transformation. It is the 
case when 9 = {0, n}, see dashed and dotted-dashed lines in figure [2] 

(ii) The second case is when there is such t% that I\{t\ + r) = —I\(t\ — r). This is the 
case of the antisymmetric driving realized for the phases 9 = {tt/2, 37r/2}, see solid and 
dotted lines in figure [2j 

The analysis of a similar system biased by a dc current and driven by an unbiased 
harmonic ac signal has been recently presented in Ref . |3T] ■ The main difference between 
these two set-ups is a constant dc current which in turn may lead to the phenomenon of 
the negative mobility (resistance) pj. Here, there is no such constant current applied to 
the system. However, by applying asymmetric ac signals we hope to be able to control 
the behavior of the first and second junctions. We would like to convince reader that 
transport across the junction can be qualitatively controlled just by adjusting the shape 
of the external signal applied to one junction only. 



4. Voltage response to biharmonic current 

To reduce a number of parameters of the model, we consider a special case of two 
identical junctions. In such a case I c \ = I c2 = I, R± = R2,a(3 = 1 and equations (pi) 
take the form 

01 = Ji(t) — sin 0i — a sin 02 + v^D £'(r), 

2 = -sin0 2 + a[/i(r) -sin0i] + VD C(r), (8) 

where the driving current I\{t) = Ji(r; a\, E\, 9) is given by the relation ([?]). 

Now, we make four general conclusions about the dc voltages as a consequence of 
the symmetry properites of the ac driving current I\{t): 

(A) - Let us consider the voltages V\ = Vi(ei) and t>2 = ^2(^1) as functions of the 
amplitude E\ of the second harmonic. If we make the transformation E\ — > —E\ to 
equations ([8]), then it follows that 

= -«i(ei), v 2 (-£i) = -v 2 (e 1 ). (9) 
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These relations yield i>i(0) = — t>i(0), 1^(0) = — 1*2(0) and we conclude that t>i(0) = 
and v 2 (0) = when the second harmonic is zero, i.e. when E\ = 0. 

(B) - The sign of voltages v\ = V\{9) and v 2 = v 2 {9) can be controlled by the phase 9. 
Indeed, if one changes the phase 9 — > 9 ± w then E\ cos(2wr + 9 ± 71) — > —s\ cos(2wr + 9) 
and from the relations (|9| one gets 

v 1 (e±n) = -v 1 (0), v 2 (9±n) = -v 2 (6). (10) 

(C) - Because the driving Ji(t) is the same for 9 — > 71 — 9 and for 9 — > 71 + 9, the 
symmetry relations 

Vl (ir - 9) = v x (ir + 6), v 2 (n - 9) = v 2 (tt + 9) (11) 

hold for any value of the phase 9. 

(D) - For the antisymmetric ac current, equations ([§]) are invariant under the time 
reversal transformation t — > —t. As a consequence [HI |2T| [28] 

Ui(tt/2) = Ui(3tt/2) = 0, v 2 (n/2) = v 2 (3ti/2) = 0. (12) 

The above set of equations ^ cannot be handled by standard analytical methods 
used in ordinary differential equations. Therefore we have carried out comprehensive 
numerical simulations. We have employed Stochastic Runge-Kutta algorithm of the 
2 nd order with the time step of [10~ 3 -j- 10~ 4 ](27r/u;). We have chosen initial phases 
0i (0) and 02(0) equally distributed over one period [0, 2ir\. Averaging was performed 
over 10 3 — 10 6 different realizations and over one period of the external driving 2ti/u. 
All numerical calculations have been performed using CUDA environment on desktop 
computing processor NVIDIA GeForce GTX 285. This gave us a possibility to speed 
the numerical calculations up to few hundreds times more than on typical modern CPUs 
[32]. 

We begin the analysis of transport properties of the system ^ by some general 
comments about its long-time behavior. In the long time limit, the averaged voltages 

(0i(r)) can be presented in the form of a series of all possible harmonics, namely, 

00 

^lim (0i(r)) = Vi + Vi(nuT), i = 1,2, (13) 

n=l 

where Vi is a dc (time- independent) component and Vi(nut) are time-periodic functions 
of zero average over a basic period. 

When both a\ and E\ are small, the dc components v 1 and v 2 of averaged voltages are 
zero in the long time limit. It can be inferred from the structure of equations pj: at least 
one of the amplitudes of the driving current I\ (r) should be greater than the amplitude 
I c i = 1 of the Josephson suppercurrent. Also for high frequency u, the averaged dc 
voltages are zero: very fast positive and negative changes of the driving current cannot 
induce the dc voltage and only multi-harmonic components of the voltages can survive. 
For smaller values of the frequency and higher values of amplitude, one can observe a 
stripe-like structure of non-zero values of both voltages v 1 and v 2 . 

We can identify four remarkable and distinct parameter regimes where: 
(I) Vx > and v 2 > 0, 
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(II) v i < and v 2 < 0, 

(III) Vi < and v 2 > 0, 

(IV) v 1 > and v 2 < 0. 

The regimes (I) and (II) dominate in the parameter space. If the regime (I) is detected 



then the regime (II) can be obtained from the relations (10) by changing the relative 
phase 9 of the biharmonic current Ii(t). Likewise, if the regime (III) is found then the 



regime (IV) can be determined from the relations (10). 
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Figure 3. (color online). Transport properties of the system of two Josephson 
junctions driven by the ac current acting on the first junction. In panels a and b 
averaged dc voltages V\ and v 2 across the first and second junction are shown for the 
coupling constant a — 0.25, temperature D — 0.001, relative amplitude of the second 
harmonic e\ = 0.5 and the relative phase 9 = 0. 



4-1. Dominated regime: the same sign of dc voltages V\ and v 2 

In figure [3] we present the regime (I) in the parameter plane {ai,tu} which illustrate 
how voltages v% and v 2 depend on the amplitude d\ and frequency u of the biharmonic 
current Ii(r) when other parameters are fixed. In this regime, both voltages always take 
non-negative values. 

The dependence of voltages V\ and v 2 on the current amplitude ai is depicted in 
panels (a) and (b) of figure |ij We observe that for smaller amplitudes of external ac 
current the transport is activated only on the first junction while the voltage across the 
second junction is zero. We can identify non-zero response of the first junction at the 
amplitude around a± — 1. The second junction awakes for larger values of amplitudes 
(the threshold is a little bit less than a\ =4). One can note non-monotonic behavior of 
all presented curves. After initial burst the amplitude of the voltage of the first junction 
decreases, revealing next enhancement together with the appearance of the non-zero 
voltage excited on the second junction (green dashed lines). For larger amplitude of the 
ac driving, voltages of both junctions decrease. In panels (a) and (b) of this figure, the 
phase 9 of the second harmonic of the external signal is different: 9 = in panel (a) 
and 9 = 71/ 3 in panel (b). The inspection of the results reveals that the change of the 
phase from 9 = 0to9 = 7i/3 reduces the dc voltages more than twice. 
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Figure 4. (color online) . The stationary averaged dc voltages v\ and v-i across the first 
(blue solid line) and second (green dashed line) junction, respectively. The parameters 
are: the dimensionless temperature D = 0.01, the frequency u> = 0.01 and strength of 
the second harmonic e\ — 0.2. The dependence on the external ac current amplitude 
a\ is presented in panels: (a) for the relative phase 9 = and (b) for 9 = 7f/3, and fixed 
coupling constant a = 0.25. In panel (c) the dependence of voltages on the coupling 
constant a is depicted for fixed amplitude ax = 5.15 and phase 9 = 0. In panel (d) the 
role of the driving current symmetry controlled by the phase 9 is shown for a = 5.15 
and a = 0.25. 



In panel (c) of figure 4, the role of coupling is illustrated: for small values of the 
coupling constant a the voltage v% is zero. There is a threshold value a c of coupling 
above which the voltage across the second junction can be generated. In the regime 
presented in figure 4, the threshold value a c ~ 0.152. From equation ^ it follows that 
the external shunt resistance R3 can be a good control parameter which decides whether 
v 2 takes non-zero or zero values. For fixed remaining system parameters, one can induce 
non-zero voltage across the second junction by decreasing the resistance R3 to values 
i? 3 < R c = (a' 1 — l)i?2- F° r R3 > R c the voltage v 2 diminishes. 

In panel (d) of figure 4, the influence of the relative phase 9 of two harmonic signals 
on the dc voltages is depicted. It is a numerical evidence for our consideration on 
symmetry relations rt9j)-(12). We can deduce two conclusions: 



(i) The maximal absolute values of voltages v\ and v 2 are generated by the symmetric 
ac current, i.e. when 9 = or 7r; 

(ii) Both dc voltages v\ and v 2 are zero when the ac current is antisymmetric, i.e. when 
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9 = tt/2 or 3tt/2. 

From Refs [T4"t [28] it follows that for small amplitudes of both harmonics of the driving 
current ([?]) the dc voltages depend on the phase as Vi oc cos# (i = 1,2). Such a behavior 
is universal and does not depend on details of systems or models. E.g. it has been 
observed in transport of flux quanta in superconducting films, see panel (b) in figure 1 
ofRef. (33]. 

As follows from equations (10), the regime (II) can be obtained from the regime 
(I) by change of the phase 9 — > 9 ± tt. Then the corresponding figures can be obtained 
from figures 3 and 4 by the reflection with respect to the horizontal axis. Therefore we 
do not discuss this regime and we do not show corresponding figures. 



4-2. Regime of the opposite sign of dc voltages V\ and V2 

Now, we address the question of whether we can identify parameter regimes (III) and 
(IV) where the dc voltages v i and v-i take opposite sign. An illustrative example of such 
a regime is depicted in figure 5. We can identify a fine stripe-like structure of regions of 
non-zero dc voltages sparsely distributed in the parameter space. Outside the stripes, 
in large regions of parameter space, the voltages V\ and Vi are negligible small or zero. 

The inspection of figures 5 suggests to perform a more accurate search for the 
horizontal section a\ = const, and the vertical section u = const.. In panels (a) and (b) 
of figure [6] we present results for the section di = 3. One can note the unique feature of 
the occurrence of windows of the frequency u where the voltage v 2 takes the opposite 
sign to the voltage v\. The voltages as a function of 00 display the non-monotonic 
dependence exhibiting maxima and minima. We can detect the interval where absolute 
values of voltages operate synchronously: simultaneous increase and decrease both of 
them. For higher temperature D = 10~ 2 and within the presented range of the base 
frequency u the voltages across both junctions alternate showing very interesting and 
rare curves of opposite voltage reversals |34J. However, if we consider ten times lower 
temperature, i.e. D = 10~ 3 , then only one window of the states working contrariwise 
survives. It means that other windows are created by thermal equilibrium fluctuations. 

Two panels (c) and (d) in figure [6] illustrate the section to = 0.1139 and the section 
uj = 0.1288 to show the influence of the amplitude a\ of the ac driving. The voltages as 
a function of a\ are non-monotonic and local maxima and minima can be detected. 

The influence of the relative amplitude E\ of the second harmonic is depicted in 
figure [7] for two values of the driving frequency to = 0.1139 (a) and u = 0.1288 (b). 
In the vicinity of Ei = 0.1 we can notice previous findings where the voltages have the 
opposite sign for both junctions. However, if we increase E\ then both voltages diminish 
and next take negative values following by positive signs for large e±, see panel (a) for 
details. For faster current u = 0.1288 the situation is opposite - for larger E\ voltages 
become positive, next negative and end up with locked state for relatively large S\. It 
is worth to mention that E\ = relates to the situation where the second harmonic is 
zero resulting in a pure symmetric driving with no possible net transport. On the other 
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Figure 5. (color online). Transport properties of the system of two Josephson 
junctions driven by the ac current acting on the first junction. In panels a and b 
dc voltages v\ and v-i are shown for the parameters set a = 0.25, D = 0.01, S\ = 0.1 
and 9 — it/ 2 revealing the most interesting response of junctions with characteristic 
stripe-like structure of regions of non-zero dc voltages. 



hand, if E\ — 1 then both harmonics of the signal ([7]) are equally strong and we are not 
able to say which part of it influences the dynamics of junctions more effectively. 

In panel (c) we show the influence of interaction between two junctions on the 
voltage characteristics. The first observation is that there are optimal values of the 
coupling strength a for which absolute values of voltages are locally maximal. The 
second observation is that variation of the coupling strength can lead to the phenomenon 
of the voltage reversal [31]: the voltage changes its values from positive (negative) values 
to negative (positive) ones. The third observation is the same as presented in panel (c) 
of figure 4: there is a threshold value of coupling a = a c (or the resistance R3 = R c ) 
such that for a < a c (or for R3 > R c ) the voltage t>2 is reduced to zero. 

In panel (d) we depict the influence of temperature. Most remarkably, transport 
is solely induced by thermal noise. Indeed, in the deterministic case (i.e. for D = 0), 
when no thermal fluctuations act, the dc voltages vanish. With increasing temperature, 
the voltage v\ starts to increase to positive- valued local maximum. Next, it decreases 
crossing zero and reaches a negative-valued local minimum. For further increase of 
temperature, the voltage V\ tends to zero. In turn, the voltage V2 takes negative values 
reaching minimal value as temperature start to increase. Next, it increases crossing 
zero and reaches local maximum. Finally, it monotonically decreases toward zero. We 
observe that optimal temperature occurs at which the voltage v\ is maximal. There is 
also another optimal temperature at which in turn the absolute value of v<i is maximal. 
Moreover, one can identify the voltage reversal phenomenon: the both voltages change 
their sign as temperature is varied. 

From the presented analysis it follows that a biharmonic drive allows one 
to conveniently manipulate transport with high efficiency by changing the system 
parameters. 
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Figure 6. (color online). The stationary averaged dc voltages V\ and V2 across 
the first (blue solid line) and second (green dashed line) junction, respectively. The 
dependence on the frequency u) of the ac driving is presented in panels (a) and (b) 
at a fixed amplitude a\ — 3.0. The temperature read D = 10 -3 on panel (a) and 
10~ 2 on panel (b). In two lower panels (c) and (d) the dependence on the external ac 
current amplitude a\ is presented at fixed frequency uj — 0.1139 (c) and to = 0.1288 (d) 
for temperature D = 10~ 2 . Parameters used in all panels are: the coupling constant 
a = 0.5, the relative amplitude of the second harmonic e± = 0.1 and the relative phase 
9 = tt/2. 



5. Summary 

This paper presents the detailed study of the system of two coupled noisy Josephson 
junctions which undergoes the influence of the external ac biharmonic current applied 
to one of the junctions only. The dependence of voltages across the first and second 
junctions upon ac driving exhibits a rich diversity of transport characteristics. In 
particular, it is possible to control both junctions to operate simultaneously with positive 
or negative voltages of the same sign or, more interestingly, with the opposite sign. 
Moreover, the phenomenon of the voltage reversals is revealed: the voltages change the 
sign as one of the parameters is varied. Our findings can be experimentally verified in an 
accessible set-up of two Josephson junctions coupled by an external resistance. The work 
can open the perspective of a new type of electronic elements which is tunable between 
positive and negative dc voltages via external control parameters like the amplitudes and 
relative phases of an ac current. The straightforward extension of the ideas presented 
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Figure 7. (color online). The stationary averaged voltages v\ and V2 across the first 
(blue solid line) and second (green dashed line) junction, respectively. The dependence 
on the relative amplitude si of the second harmonic of the ac driving is presented in 
panels (a) and (b) for fixed coupling a = 0.5 and for two frequencies w = 0.1139 (a) 
and uj — 0.1288 (b). In panel (c) the dependence of voltages on the coupling constant 
a is depicted for e% — 0.2 and 10 = 0.12884. In panel (d) the role of temperature 
controlled by the dimensionless parameter D is shown for E\ = 0.2, u; = 0.1288 and 
a = 0.5. Parameters used in all panels read: relative phase 9 = it/ 2 and amplitude 
ai = 3.0. In panels (a)-(c) temperature D = 10 -2 . 



here would cover an increasing number of coupled junctions to three (or more) as well 
as study of other types of driving. 
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